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Abstract 

Consider a layer made of a m3m insulator crystal, with faces cut 
parallel to a symmetry plane. Then bond it onto a semi-infinite mm2 
piezoelectric substrate. For a X- or Y-cni of the substrate, a Love 
wave can propagate in the resulting structure and the corresponding 
dispersion equation is derived analytically. It turns out that a fully 
explicit treatment can also be conducted in the case of a y-cut rotated 
about Z. In the case of a germanium layer over a potassium niobate 
substrate, the wave exists at any wavelength for X- and y-cuts but 
this ceases to be the case for rotated cuts, with the appearance of 
forbidden ranges. By playing on the cut angle, the Love wave can be 
made to travel faster than, or slower than, or at the same speed as, 
the shear bulk wave of the layer. A by-product of the analysis is the 
derivation of the explicit secular equation for the Bleustein-Gulyaev 
wave in the substrate alone, which corresponds to an asymptotic be- 
havior of the Love wave. 

1 Introduction 

Layered structures, especially film/coating substrate systems, play an im- 
portant role in micro-electro-mechanical systems (MEMS) and in microelec- 
tronics packages. In order to achieve high performance, many surface acous- 
tic wave (SAW) devices/sensors are made of layered architectures such as, 
for instance, a dielectric, or a piezoelectric, or a non-piezoelectric semicon- 
ductor (finite-thickness) layer deposited onto a (semi-infinite) piezoelectric 
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substrate. For certain configurations, it is possible to have a one-component 
wave travel in the structure, in the direction of the interface: this guided 
(shear-horizontal) Love wave leaves the upper face of the layer free of me- 
chanical tractions, its amplitude varies sinusoidally through the thickness of 
the layer and then decays rapidly with depth in the substrate, and it is such 
that all fields are continuous at the layer /substrate interface. 

Love waves in piezoelectric layered acoustic devices are most suitable for 
high frequency filters because of their high phase velocity, and they also show 
great promise in bio-sensors applications with liquid environments because 
of their high sensitivity. Consequently, they have received much attention 
over the years. For example, Lardat et al. [1] found under which conditions 
a piezoelectrically stiffened Love wave exist and presented experimental and 
analytical results on surface wave delay lines. Kessenikh et al. [2] investigated 
surface Love waves in piezoelectric substrates of classes 6, 4, 6mm, 4mm, 622, 
and 422 with an isotropic dielectric layer. Hanhua and Xingjiao [3] studied 
Love waves for a structure made of a 6mm piezoelectric layer and a 6mm 
piezoelectric substrate, with a common symmetry axis in the plane of the 
interface, while Darinskii and Weihnacht [4] had a similar structure, made 
of 2mm piezoelectric layer and substrate with common symmetry axes, one 
of which is aligned with the propagation direction and another is aligned 
with the normal to the interface. Jakoby and Vellekoop [5] reviewed the 
properties of Love waves and associated numerical methods for a piezoelec- 
tric/piezoelectric layered structure, when the substrate is made of ST-cut 
quartz; so did Ogilivy [6], with special emphasis on the mass-sensitivity load- 
ing of biosensors. The reader can find additional pointers to the literature 
on Love waves for piezoelectric sensors in those articles and in the references 
therein, as well as in the reviews by Farnell and Adler [7] or by Gulyaev [8]. 

The present work is concerned with the propagation of Love waves in 
a composite structure, made of a m3m cubic insulator (non-piezoelectric 
semiconductor) layer of finite thickness /i, bonded onto a Y-cut rotated about 
Z, mm2 orthorhombic, piezoelectric substrate, see Fig{Tl With this singular 
structure, it is possible to exploit the potentiality offered by an efficient direct 
analytic method. The steps to be followed are: (i) Determine the general 
solution for the layer, satisfying the boundary conditions on the upper face; 
(ii) Derive some fundamental equations in the substrate using the Stroh 
formulation of the equations of motion there; (in) Match the solutions at 
the layer /substrate interface; and (iv) Deduce the dispersion equation for 
piezoelectric Love waves in explicit form. 
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The dispersion relation admits a denumerable number of solutions (bran- 
ches) and the solutions for the mechanical displacements, shear stresses, elec- 
tric potentials, and electric inductions are deduced explicitly in the case of 
a metalized, mechanically free, upper surface, brought to the zero electric 
potential (short-circuit). In passing, the explicit secular equations for the 
Bleustein-Gulyaev wave speed and for the limiting wave speed (substrate 
only, no overlayer) are found as a cubic and as a sextic in the squared wave 
speed respectively, for rotated cuts. The special cases of a X-cut or F-cut are 
treated separately. The effects of the angle of the cut on the phase velocity 
of the first modes are illustrated numerically for a specific layered structure, 
namely a germanium layer over a potassium niobatc substrate and the ap- 
pearance of a forbidden band of frequency is uncovered for a rotated cut, 
in sharp contrast with the non-rotated cuts where the waves exist for all 
frequencies. 
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2 The layer 

First consider the upper layer, which is made of a cubic, m3m non-piezoelectric 
(semi-conductor) crystal, with mass density p. For two-dimensional motions 
(independent of xs), the anti-plane equation of motion decouples from its 
in-plane counterpart and reads (0 < 0:2 < h) 

C44(m3,11 + ^3,22) = PM3,«, (1) 

where Ms is the anti-plane mechanical displacement, C44 the transverse stiff- 
ness, and the comma denotes partial differentiation. For a solution in the 
form of an inhomogeneous wave traveling with speed v and wave number k 
in the xi direction, such as 

^3 = f/3(A;x2)e''=(^^-^*), (2) 
(where t/3 is a function of kx2 alone), the equation of motion ([!]) reduces to 

f>3 + - 1) = 0, where v := ^/t^p. (3) 
The general solution to this second-order differential equation is either (i) 



Uz{kx2) = f/3(0) cos \/ - Ikx 



y2 



^2 



y2 



y2 



-1-^4 sin \l — — lkx2 ] , when v > v, (4) 

or (ii) 



Us{kx2) = UsiO) I cosh ^/l - ^kX2 



-|-y4sinh \/ 1 — ^kx2 | , when v < v, (5) 



where A is a constant to be determined from the boundary condition at 
X2 = h. This latter condition is that the upper face of the layer be free 
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of mechanical tractions, so that (T23 = £44^3^2 is zero there. Then it follows 
from (jl]) and ([5]) that: A = tan a/ f ^/ v'^ — Ikh in Case (i), and that: A = 
— tanh \/l — v'^/v^kh in Case (ii). Consequently the mechanical field is now 
entirely known in the layer. 

The procedure to find the electrical field is similar, and even simpler 
because the layer is not piezoelectric. Taking the electric potential in the 
form 

= if{kx2)e''^''-'"\ (6) 
(where is a function of kx2 alone), Poisson's equation: Acf) = reduces to 

0"-0 = 0, (7) 

with general solution 

ip{kx2) = (^(0) (cosh kx2 + B sinh kx2) , (8) 

where i? is a constant to be determined from the boundary condition on the 
upper face of the layer X2 = h. A thin metallic film covers that face, and 
it is grounded to potential zero: 0{kh) = 0. It then follows from ([H]) that 
B = — coth kh, and consequently, that the electrical field is now entirely 
known in the layer. 

For the problem at hand, only the values of the fields at the interface X2 = 
between the layer and the substrate are needed. Because the mechanical 
displacement U3 and the electrical potential are in the forms ([2]) and ([6]), 
respectively, the mechanical traction (T23 and the electrical displacement D2 
must also be of a similar form, due to the constitutive relations: 6"23 = €44113^2 
and D2 = — eii0,2, where en is the dielectric constant of the layer. Hence, 

^23 = iki23{kx2)e"'^''-'''\ D2 = iA;4(A;x2)e*^(^^-''*), (9) 

(say), where £23 and d2 are functions of kx2 alone. In particular, the conclu- 
sions drawn from the calculations conducted above are that at the layer/sub- 
strate interface, 

t23(0) = -icUsiO), 4(0) = -ieV(O), (10) 
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where 



c := —C44\ 1 — TTT tanh \ 1 — ^kh, when v < v, 



c :— C44y — — Itan y ^ — ikh, when v > v, 




e :— ill cothkh. 



(11) 



3 The substrate 



The substrate occupies the half-space 0:2 < and is made of the superstrong 
piezoelectric crystal, potassium niobate KbNOa, with orthorhombic mm2 
symmetry. Cut the crystal along a plane containing the Z axis and making 
an angle 9 with the XY plane. Let X1X2X3 be the coordinate system obtained 
after the rotation 



Yet again, for a two-dimensional motion (independent of X3) the anti-plane 
strain and stress decouple from their in-plane counterparts. Hence, with 

Ui = U2 = 0, Ma = u^{xi^X2^t) for the mechanical displacement, and = 
(j){xi,X2,t) for the electric potential, the constitutive relations yield an = 
C22 = CTas = C12 = for the stress components and = for the electrical 
displacement, and they reduce to 



where the Cij, Cij, 6^ are related to the corresponding quantities Cij, Cij, iij in 
the crystallographic coordinate system XYZ through tensor transformations 



m n 
~n m 
1 



where m — cos6, n — sm6. 



(12) 



D2 



C44U3,2 + C45U3,1 + 6140,1 + 6240,2, 
C45W'3,2 + C55M3,1 + 615^,1 + ei40,2, 
614^3,2 + 615^3,1 - 6110,1 " 6120,2, 
624^3,2 + 6i4M3,l " 6l20,l - 6220,2, 



(13) 
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[9] as: 



C44 = rn^cu + n^c^^, 



C45 = mn[cu - C55 
ei5 = 
en = 
ei2 = mn(e22 - hi 



C55 = n^C44 + m^C55, 
624 = '"^^624 + ^^ei5, 

ei4 = mn(e24 - 615), 
€22 = ri^eii + m^e22, 



(14) 



Now consider the a:2-cut, Xi-propagation of a Shear-Horizontal interface 
acoustic wave that is, a motion with speed v and wave number k where the 
displacement field M3 and the electric potential are of the form, 



{M3,0}(xi,X2,t) = {t/3(fcx2),</^(A;x2)}e^'^("^-^*), 



(say), with 



U2.{-oo) = 0, v'(-oo) 



0. 



(15) 
(16) 



It follows from the constitutive equations (fT3l) that the traction (732 and the 
electric induction D2 are of a similar form. 



W32, D2}{xi,X2, t) = ik{ts2{kx2), d2{kx2)}e 



ik{xi—vt) 



(say), where 



(17) 
(18) 

The non-trivial part of the equations of piezoacoustics, aijj = pui^tt, 
Di,i = (where p is the mass density of the crystal), can be written as a 
second-order differential system [10]: 



t32(-00) = 0, rf2(-00) 



0. 



T 



+ 2iR 



-[Q-{pv')J] 



'Us 




"0" 










(19) 



where the prime denotes differentiation with respect to kx2 and 

rp ._ C44 624 ^ 

[e24 —£22] ' 

Q := h 1 ^ J ■■- 

ei5 -en 



C45 ei4 
ei4 -ei2 
1 0" 







or as a first-order differential system in the form 

^' = iA^I, where ^ikx2) := [f/3, ^31, ^^2]*, 



(20) 



(21) 
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and has the Stroh block structure 



N 



Ni N2 



with [11] 



-T-^R, N2 = T-\ N3 = RT-^R - Q. 



(22) 



(23) 



Here the components of are easily computed from (120|) and (l23!l . by hand 
or using a Computer Algebra System. 

Seeking a solution to (HM in the form: [U3, cpY = [Us{0), ip{0)Ye^'"^^'^, where 
q is constant, yields the propagation condition 



det [q^T + 2qR + Q- {pv^)J] = 0, 



(24) 



a quartic in q. Explicitly, 

+ 2^4g^ + {eg - pv')q' 



^22 



+ 2^(4-p.^)g+^(cS-p.^) 

^22 ^22 



0, (25) 



where 



-44 



-45 



-55 



C44 + 



^24 D ,^22 „ 615624 

— , qg = C44 + C45 h 2 , 

622 612 ^12 

, en 615624 612 , .614 
C55 + C44 h 2 h 4645 h 4 — , 

622 622 622 622 



655 + 



-15 
611' 



-26 



I 611 „ 615624 

655 + 644 h Z . 

612 612 



(26) 



Out of the four possible roots, only two have a negative imaginary part, in- 
suring exponential decay as X2 — > —00. Computing these qualifying roots 
analytically is not an easy matter, because here the speed v is still an un- 
known. Although it is actually possible to do so [12,13], this paper follows 
a different route and makes extensive use of the "fundamental equations" 
derived in [14-16]. They read 



^(0) ■ M(")|(0) = 0, where M^") 
and n is any positive or negative integer. 



/ 
/ 



AT", 



(27) 
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4 The layer/substrate structure 



The continuity of all fields at the layer/substrate interface X2 = imposes 
the boundary conditions: 

^23(0) =£23(0), d2(0)=4(0). (28) 

Now the dispersion equations are derived explicitly. In the special cases of X 
and F-cuts, the dispersion equation is exact: if it is satisfied, then the Love 
wave exist. In the other cases, the dispersion equation is rationalized: it has 
spurious roots, not corresponding to a true Love wave so that a subsequent 
validity check is necessary. 



4.1 Special cases ^ = 0, 90° 

For a X-cut or a Y-cut of the substrate, the analysis simplifies considerably 
and a direct treatment is possible, leading to an exact dispersion equation, 
and not requiring the use of the fundamental equations. 

When 771 = 1,71 = 0, orm = 0,?7, = l, the parameters C45, e^, 612 vanish 
according to (fl^ . Then in fl20|) . R = also, and the quartic (p4l) becomes 
the following biquadratic [4]: 

q^-Sq^ + P = 0, (29) 
where the non-dimensional quantities S and P are given by 

cu^n + 055622 + 2615624 - e22{pv'^) 



-S 



644622 + 634 



P = '^^'n + ^l^-j(P-')_ (30) 
644622 + 624 

The relevant roots qi and q2 come in one of the two following forms, either 
(a) qi = -i/?i, q2 = -i/?2, 

or (6) qi = —a — i/?, q2 = a — i/3, (31) 

where > 0, /5 > 0. In either case, qi + q2 has no real part and a negative 
imaginary part, and qiq2 is a negative real number. Explicitly, 



qi + q2 = -iV - S, qiq2 = -VP. (32) 
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The associated eigenvectors <^"'" and C,"^ follow from (for instance) the third 
column of the adjoints to the matrices — qil and — g2-^? as = [a*, 6*]*, 
(i = 1, 2) where 



a 



2 , ^24 2 , ^15 



^22 £22 



^22 



(33) 



(here T is given by ( 120|) . with components evaluated at 0° or 90°). Then the 
general solution to the equations of motion is of the form 



|(A;x2) = 7ie''''^"^C' + 72e''''^"^C^ 



(34) 



where 71, 72 are disposable constants. 
At X2 = 0, it can be split as 



[[/3(0), v^(O)]* = A7, [t23(0), d,m' = B-f, 



where 



A:=[a^\a% B := [b'\b% 7:= [71,72]*. 
Now the boundary conditions fl28|) and ffTOj) give the link 

57 = [^23(0), 4(0)]* = Diag(-ic,-ie)[f/3(0),<^(0)]* 
= -i Diag(c, e)A7, 

from which the dispersion equation follows as 

|iS^-^ -Diag(c,e)| = 0. 



(35) 
(36) 



(37) 



(38) 



It is written in this form to take advantage of the many desirable properties 
of the surface impedance tensor iBA~^ [24]; thus, this matrix is Hermitian 
in the subsonic range (defined below), and of the compact form: 



iBA 



-1 



qi + q2 



PV - C55 - C44gig2 -ei5 + 62451^2 



-ei5 + e24gig2 



cii — ^22qiq2 



(39) 



Moreover, the eigenvalues of the aggregate impedance tensor in (13811 are real 
monotonic decreasing functions of v for any fixed kh, so that the wave speed 
of each mode is obtained unambiguously from the roots of (l38l) . see Shuvalov 
and Every [25]. 
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Using the identities qf + q2 = S, = P and the connections (13U]1 and 
the exact, explicit, dispersion equation is finally derived as 



^ - C55 - , . ei5 + 624^^ 

/ ^ ^ c . 

ei5 + e24V^ en + 622 , , 

-U /: 



2/P- 5 



2/P-5 



(40) 



which is a fully explicit equation, because c, e are defined in ffTOj) and 5, P 
in ([30]). 

The dispersion equation is valid in the subsonic range that is, as long 
as the speed is less than the limiting speed vl, the smallest speed at which 
the biquadratic fl2^ ceases to have two roots gi, q2 with a positive imaginary 
part. When qi, q2 are of the form (a) m flM]) . then is root to: P = 0; when 
they are of the form (b), then vl is root to: 2\/P = S. In either case, the wave 
becomes homogeneous at v = vl because then the roots of the biquadratic 
are real; they are: ±g, where q = y/Si or q = a/ Sl/2, according to each 
case (here Sl is S given by (|5UI) at v = vl). The associated eigenvectors are 
[a, ±6]* where 



a 



2 , en 624 2 I <2l5 

g + — , — q + — 



£22 £22 



^22 



b = qTa. 



(41) 



The vanishing of the wave away from the interface can no longer be insured 
then, but the continuity of the fields at X2 = can. The conclusion is that 
the boundary conditions (!28l) and (ITOl) lead to the following cut-off equation, 



Caa + 



"24 
^22 



2 , £44611 + 624615 

g H — 



£22 



en 

624 ei5 

£22 



2 I ^11 

q^ ^ 

^22 



624 2 , ^15 

— q^ H 

^22 £22 



(42) 



This equation has an infinity of roots in kh (corresponding to the intersections 
of the graph of tan with the graph of coth). Each root {kh)L say, is a cut-off 
parameter for each dispersion mode, at which the Love wave ceases to exist. 

Note that the dispersion equation fHOj) is consistent with the secular equa- 
tion of a Bleustein-Gulyaev wave traveling in the substrate alone: as h tends 
to zero, c — >• 0, 6 — > 00 so that it reduces to: pv"^ — C55 = C44A/P, which, once 



11 



squared, coincides with the quadratic (15^ . obtained in the next Section. It 
is also consistent with the dispersion equation for a purely elastic Love wave. 
Indeed, by taking Cij and eij — > in fl30|) and fHOj) . the equation of 
Lardat et al. [1] is recovered: 




C44 
C44 



C55 - pv"^ 



tan \l^^—^kh= ' I . (43) 

/ pv^ — C44 

C44 



Finally, it is consistent with the dispersion equation of Love surface waves 
in an isotropic dielectric layer over a 6mm piezoelectric substrate [2] , by the 
corresponding specialization. 



4.2 Rotated cut 

Combining the boundary conditions (128!) with the results for ^23(0) and (^2(0) 
of ( ITOl) gives the following form for ^(0): 

^(0) = f/3(0)[l,a,-ic,-iea]*, (44) 

where a := ip{0)/U3{0) is complex: a = ai + ia2, say. Then the fundamental 
equations (p7l) read 



1 " 








Mis' 






1 




a 














a 


= 


ic 










JW34 




— ic 


iia 






JW24 




JW44 _ 




—iia 





(45) 



or 

[Mg^ + ceM(?](2ai) + [eM^^^ - cM(3")](2a2) 

+ [Mi? + e2Ml?](a? + a^) = -[mJ? + c'm!^^ (46) 

Writing them for n = —1, 1, 2, and re-arranging the three resulting equations, 
leads to the following non-homogeneous system of linear equations, 

[fcl|fc2|fc3]P= -fc4, (47) 



12 



where p := [2,ai,2a2,Oil + ^2]* fci, fc2, ^3, ^4 are the vectors with 
components: 



Mg^+e^Mi:^ m};^) + c2Mi3"\ (48) 



(n = —1,1,2) respectively. Cramer's rule gives the unique solution to the 
system as 

2ai = -Ai/A, 2a2 = -A2/A, aj + aj = -A3/A, (49) 

where A := det[A;i|fc2|fc3]j Ai := det [A;4 1 ^2 1 ^3] , A2 := det[fci 1^41^3], and 
A3 := det[fci|fc2|fc4]. The dispersion equation follows then from the compat- 
ibility of the equalities (H^ : 

Al + Al + 4A3A = 0. (50) 

When (and if) this dispersion relation yields a real positive wave speed 

V for a given wave number k, it remains to be checked whether that speed 
corresponds to a valid solution. Proceed as follows. First recall that the 
exact boundary condition is of the form fl38l) . where now the a*, 6* {i = 1, 2) 
are defined by 

i 2 I 0^12 , en 624 2 I 0*^14 , ei5 

€22 ^22 ^22 ^22 ^22. 

:= {q{r + R)a\ (51) 

The computation of the corresponding surface impedance tensor iBA~^ is 
long but perfectly possible analytically; its components depend on gi, q2 
through the sum qi + q2 and the product qiq2- Now, having found a speed 
from (!50!) . compute numerically the roots of the quartic fl2^ . Select qi and 
g2, the roots with negative imaginary parts (if there are no such roots, then 

V is not valid.) Then compute qi + g2, 9192, and iBA~^, and check whether 
the exact boundary condition (1551) is satisfied. 

Finally, it is also possible to determine exactly the limiting speed at 
which the imaginary part of an attenuation factor first vanishes. Fu [12] 
shows that vl is the smallest root of D = 0, where D is the discriminant of 
the cubic resolvent associated with the quartic ( 12^ . Thus rewrite the quartic 
( 12^ in its canonical form, 

+ rp^ + sp + t = 0, (52) 
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say, using the substitution q = p — (l/2)(ei2/e22)(cf^/c^)- Then, 

D = -4(12t + + (-72tr + 2r^ + 275^)^ (53) 

Here the equation D = turns out to be a sextic in the squared wave speed, 
according to Maple. 



5 Bleustein-Gulyaev wave as h ^ 

Shuvalov and Every [25] show that a great variety of asymptotic behaviors 
arises for an interface wave on a coated half-space, depending on whether the 
layer or the substrate is "fast" or "slow" , or "dense" or "light" , depending on 
which mode is under consideration, and depending on several other factors. 

When the thickness of the layer vanishes here, the asymptotic behav- 
ior of the Love wave in the layer /substrate structure is that of a shear- 
horizontal surface wave propagation over the piezoelectric substrate alone 
(the Bleustein-Gulyaev wave [17, 18]). For such a wave, with metalized 
boundary conditions, the vector ^(0) takes the form 



m = PsiO), 0, 0, c?3(0)]* = f/3(0)[l, 0, 0, a] 



(54) 



where a := (i3(0)/[/3(0) is complex. The fundamental equations (127|) . written 
for n = —1, 1, 2 (say), can be arranged as 



r{2) 



M- 



14 



M, 



(2) 



44 





1 




"0" 




a + a 









aa 








(55) 



a homogeneous linear system of three equations. Its solution is non-trivial 
when the determinant of the 3x3 matrix on the left hand-side is zero. The 
resulting secular equation is a cubic in pv"^. At ^ = and 9 = 90°, the 
secular equation factorizes into the product of a term linear in pv^ and a 
term quadratic in pv^. In particular, at ^ = the quadratic is: 



[pV -C55) +C44 



C55jeil 



^15 



044622 + ^24 



0. 



(56) 



Now for KbN03, the material parameters of interest are [19]: C44 = 7.43, 
C55 = 2.5 (10^0 N/m2), g24 = 11.7, ^15 = 5.16 (C/m^), e~n = 34eo, €22 = 780eo 
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(eo = 8.854 x 10 F/m), and p = 4630 kg/m^. Using these values, the 
corresponding parameters in the rotated coordinate system follow from fll4p . 
and in turn, T, R, Q follow from ([20]), from ([22D-([23D, and M^") from ([27]). 
Then the cubic secular equation is solved for v for any value of the cut angle 
6. Out of the three possible roots, only one may correspond to the Bleustein- 
Gulyaev wave (Ref. [20] explains how the adequate root is selected.) It turns 
out that the wave exists for all angles, with a speed vbg (say) increasing from 
2895.35 m/s at 6' = to 4450.85 m/s at 9 = 90°. Figj2] shows the dependence 
in 6, and is in agreement with the plot obtained by Nakamura and Oshiki [21] 
and by Mozhaev and Weihnacht [22]. 

FiglD also displays the speed of the bulk shear wave in a layer made 
of germanium, for which [23]: C44 = 67.1 x 10^° N/m^, in = 16.6eo, and 
p = 5330 kg/m^; here -0 = 3550.31 m/s. The angle at which vbg = "0 is 
^0 = 50.0817°. Finally, the variation of the limiting speed vl in KbNOs with 
the angle of cut is shown as well, and that plot is also in agreement with 
Mozhaev and Weihnacht [22]. 

6 Dispersion curves 
6.1 Special case 6 = 0° 

At 6 = 0° and kh = 0, the interface wave travels in the substrate alone, 
with the Bleustein-Gulyaev wave speed of 2895.35 m/s. The limiting speed 
(found here as root of = S) is vl = 3939.33 m/s. The speed of the 

fundamental mode starts at the Bleustein-Gulyaev wave speed at kh = 0, 
increases to a maximum speed of about 3857.18 m/s at kh = 0.517, and 
then decreases toward the shear bulk speed of the layer, v = 3550.31 m/s. 
In the narrow range where kh is smaller than 7.08436 xlO~^, the speed v of 
the fundamental mode wave is smaller than v and thus c is given by ([TT]) 9: 
otherwise it is given by ([TTl) i . The fundamental mode exists for all values of 
kh. 

The speeds of the subsequent modes start from the limiting speed f/, at 
{kh)L and tend to t) in a monotone decreasing manner. The cut-off parameter 
{kh)L for the first mode, second mode , and third mode is: 6.24, 12.75, and 
19.27, respectively. Figl3] shows the dispersion curves for the fundamental 
mode and for the first mode. Qualitatively, the plots echo those of Kielczynski 
et al. [26] who considered a structure made of a 6mm substrate covered with 
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Fig. 2: Solid curves: Variations of the Bleustein-Gulyaev wave and limiting wave 
speeds with cut angle in a homogeneous KbNOs substrate. Dashed line: Speed of a 
(bulk) shear wave in germanium. 
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a "depolarized" layer. 

6.2 Special case e = 90° 

At = 90° and kh = the limiting speed is root of P = 0; here vl = 4508.73 
m/s. The speed of the fundamental mode starts at the Bleustein-Gulyaev 
wave speed of 4450.85 m/s at kh = 0, increases to a maximum speed of about 
4480.69 m/s at kh = 0.0581, and then decreases toward the shear bulk speed 
of the layer, v = 3550.31 m/s. Here too, the fundamental mode exists for all 
values of kh. 

The speeds of the subsequent modes start from the hmiting speed vl at 
{kh)L and decrease toward -0. The cut-off parameter (kh) l for the first mode, 
second mode , and third mode is: 4.05, 8.06, and 12.06, respectively. Figl3] 
shows the dispersion curves for the fundamental mode and for the first and 
second modes. 

6.3 Special case = Oq 

As an example of a rotated cut, consider the case where the speed of the 
Bleustein-Gulyaev wave in the substrate is equal to the shear wave speed 
in the layer; this occurs at 6* = 6*0 = 50.08°, see Figj2l For this cut the 
limiting speed is vl = 3813.36 m/s. Starting from vbg = v at kh = 0, the 
speed of the fundamental mode increases rapidly with kh. At kh = 0.1131, 

V = Vl and the wave ceases to exists. This state of affairs continues until kh 
reaches 0.9410, after which the wave exists and its speed decreases toward 
V. Hence a forbidden band of frequencies emerges for the fundamental mode, 
in clear contrast with the situation for non-rotated cuts. Note that the 
dispersion equation fl5U]) actually gives roots below vl in that range, which 
must nevertheless be discarded as they do not satisfy the exact boundary 
condition (!38|) . 

Here the first mode starts at the cut-off parameter {kh)^ of: 8.946 with 

VI and then decreases toward v. FigjH provides a zoom into the dispersion 
curves of the fundamental mode around the forbidden band and of the first 
mode up to kh = 12. 

In this example, the layered structure supports a shear-horizontal wave 
which in the long and short wavelength ranges travels with the speed of the 
layer's bulk shear wave; in the intermediate range, the wave either does not 
exist, or travels at a greater speed. 



17 



VBG (90°) 



VL (90°) 



4400 



4200- 



4000 



3800- 




3600 



3400- 



3200 



3000- 



Fig. 3: Dispersion curves at 6* = 0°, 6*0, 90°. 
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